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. Abstract 

' We show how to search N items arranged on a x grid in time 0{\/N log N), using 

. a discrete time quantum wallc. This result for the first time exhibits a significant difference 

^ 1^ I between discrete time and continuous time walks without coin degrees of freedom, since it has 

P5 ' been shown recently that such a continuous time walk needs time D,{N) to perform the same task. 

^ , Our result furthermore improves on a previous bound for quantum local search by Aaronson 
and Ambainis. We generalize our result to 3 and more dimensions where the walk yields the 

^ ' optimal performance of 0{y/N) and give several extensions of quantum walk search algorithms 

• i-H , for general graphs. The coin-flip operation needs to be chosen judiciously: we show that another 

' "natural" choice of coin gives a walk that takes Vl{N) steps. We also show that in 2 dimensions 

■ it is sufficient to have a two-dimensional coin-space to achieve the time 0{-/N log N). 



1 Introduction 

Quantum walks are quantum counterparts of classical random walks. Classical random walks have 
many applications in randomized algorithms |MR,95j and we hope that quantum walks would have 
similar applications in quantum algorithms. Both discrete-time |Mey96[ lAAKVOTl lABN^oT] and 
continous time |FG98l CFG02 quantum walks have been introduced^. The definitions of the two 
are quite different. In continous time, one can directly define the walk on the vertices of the graph. 
In discrete time, it is necessary to introduce an extra "coin" register storing the direction in which 
the walk is moving. 

Because of this difference in the definitions, it has been open what the relation between discrete 
and continous walk is. In the classical world, the continous walk is the limit of the discrete walk 
^For an introduction to quantum walks see |Kem03a| . 



1 



when the length of the time step approaches 0. In the quantum case, this is no longer true. Even 
if we make the time-steps of the discrete walk smaller and smaller, the "coin" register remains. 
Therefore, the limit cannot be the continous walk without the "coin" register. This means that one 
variant of quantum walks could be more powerful than the other in some context, but so far all 
known examples have given similar behavior of the two walks (see e.g. |( Jb\ ^021 IKemOS bl lCCD+n3j ) . 

In this paper, we present the first example where the discrete walk (with "coin") outperforms 
the continous walk (with no "coin" ) . Our example is the spatial search |Benn2| lAAOS) variant of 
Grover's search problem. In the usual Grover's search problem |(irro96j . we have N items, one of 
which is marked. Then, we can find the marked item in 0{\/N) quantum steps, with one quantum 
step querying a superposition of items. In contrast, classically ^}{N) queries are required. In the 
"spatial search" variant, we have the extra constraint that the N items are stored in N different 
memory locations and we need time to move betwen locations. This may increase the running time 
of a quantum algorithm. 

The first "spatial" version of Grover's algorithms with optimal performance was given by 
|SKW03] who showed how to search items arranged on the n-dimensional hypercube, using 
a discrete quantum walk. 

In this paper, we consider the 2-dimensional arrangement where N memory locations are ar- 
ranged in an \/iV x ^/N grid. This was first studied by Benioff |Ben02j who observed that the usual 
Grover's search algorithm takes i}{N) steps. It uses B(\/iV) query steps but, between each two 
queries, it might move a distance of Q{VN). Thus, the total time becomes Q{N) and the quantum 
speedup disappears. Aaronson and Ambainis |AA03j fixed this problem by giving an algorithm for 
searching the 2-dimensional grid in 0(\/iVlog^ A^) total steps^ (counting both queries and moving 
steps) and the 3-dimensional grid in 0(\/iV) steps, using Grover's algorithm together with multi- 
level recursion. Quantum walks were first applied to this problem by Childs and Goldstone |CG03j 
who studied the search on the grid by a continous quantum walk. They discovered that the conti- 
nous walk provides an alternative search algorithm with optimal performance of 0{^/N) in 5 and 
more dimensions, but not in 2 or 3 dimensions, where the continuous walk takes Q{N) and Q{N^/^), 
respectively. In 4 dimensions the continuous time walk algorithm performs as 0{^/N log N). 

In this paper, we use discrete-time quantum walks to design an algorithm that searches the 
grid in 0{\/lV log N) time in 2 dimensions and 0{^fN) time in 3 and more dimensions. Thus, our 
algorithm is faster than both the non-walk quantum algorithm of |AA03j and the algorithm based 
on the continous time quantum walk |GG03j . In addition to having a very simple structure our 
algorithm also uses only 1 or 2 qubits of extra memory (or log2d qubits for the d-dimensional grid) , 
besides the current location. (The previous algorithm of |AA03j uses O(log'^n) qubits of extra 
memory.) 

Besides improving the running time, we present several interesting features of quantum walks. 
The first feature is that the discrete-time walk succeeds while the continous walk does not. Secondly, 
the behavior of the discrete quantum walk on the grid crucially depends on the choice of the "coin" 

^The running times for the 2-dimensional grid are for the case when the grid contains one marked item. The 
general case (an arbitrary number of marked items) can be reduced to the one item case with a log A'^ increase of 
the running time |AA03| . That would result in a running time of 0(y/N log^ N) for the algorithm of |AA03| and 
0{y/N \og^ N) for our algorithm which we present in this paper. For 3 and higher dimensions, the general case can 
be reduced to the one item case with just a constant factor increase |AA03| . Thus, the asymptotic running times 
stay the same. 
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transformation. One natural choice, discovered numerically by Neil Shenvi |Slie03j . leads to our 
algorithm while some other natural choices fail to produce a good algorithm. Thus, the "coin" 
transformation could be a resource which affects the algorithm profoundly. We give both upper 
and lower bounds for the performance of some natural choices of the "coin" . Surprisingly we show 
that in the case of the 2-dimensional grid only 2 (and not the standard 4) coin-degrees of freedom 
are sufficient to achieve the quantum speed-up. The insights gained from our study might aid in 
the design of future discrete quantum walk based algorithms. Several such algorithms have recently 
been discovered jCCD+OSl lAmbn.SI IMSSn.SL KJhmi ISzen4j . 

Our presentation allows a fairly general approach to quantum walk search algorithms on graphs. 
In particular it simplifies the proof of jSKWOSj , where the relevant eigenvectors had to be "guessed" . 
We also give a discrete walk search algorithm on the complete graph and show its equivalence to 
Grover's algorithm and outline several generalizations of our results. 

2 Preliminaries and Notation 

2.1 Model 

Our model is similar to the one in |AA03j . We have an undirected graph G = {V, E). Each vertex v 
stores a variable S {0, 1}. Our goal is to find a vertex v for which = 1 (assuming such vertex 
exists). We will often call such vertices marked and vertices for which Uy = unmarked. 

In one step, an algorithm can examine the current vertex or move to a neighboring vertex in 
the graph G. The goal is to find a marked vertex in as few steps as possible. 

More formally, a quantum algorithm is a sequence of unitary transformations on a Hilbert space 
Ti.i0Ti.v- Tiy is a Hilbert space spanned by states \v) corresponding to vertices of G. Tii represents 
the algorithm's internal state and can be of arbitrary fixed dimension. A i-step quantum algorithm 
is a sequence C/i, U2-, ■ ■ ■, Ut where each Ui is either a query or a local transformation. A query Ui 
consists of two transformations {U^, U}). Uf ® I is applied to all Tii® \v) for which = and 
Ul / is applied to all Tii(® \v) for which = 1. 

A local transformation can be defined in several ways |AAn3j . In this paper, we require them 
to be local. A transformation Ui is Z-local if, for any v £ V and \iJj) G Tii, the state Ui{\^p) <® \v)) 
is contained in the subspace Tii ® '^r(?)) where 7Yr(t;) C Tiy is spanned by the state \v) and the 
states \v') for all v' adjacent to v. Our results also apply if the local transformations are C-local 
(another locality definition introduced in .AAO^)- 

The algorithm starts in a fixed starting state \ip start) and applies Ui, . . ., Ut- This results in a 
final state \ip final) = UtUt~i - ■ ■ Ui\ipstart) ■ Then, we measure l^Jstart)- The algorithm succeeds if 
measuring the Tiy part of the final state gives \g) such that ag = 1. 

For more details on this model, see 

2.2 Search by quantum walk 

In what follows we will assume that G is undirected and d-regular, i.e. has constant degree d. To 
each vertex we can associate a labeling {1, . . . , d} of the d edges (directions) adjacent to it and an 
auxiliary "coin" -Hilbert space Tid = . . . , \d)}. Let Ti]\f be the Hilbert space spanned by the 
vertices of the graph, then the walk takes place in the joint space of coin and graph Ti = Tid ^3 Ti]\f. 
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Definition 1 [Discrete Quantum Walk on G:] The discrete quantum walk is an alternation of coin 
flip and moving step: U = S ■ C, where S is a shift controlled by the coin register 

S : |z) (g) \x) — \Tr{i)) (g) \x) (1) 

i = 1,. . . ,d and x,x & V, x and x are connected by the edge labelled "i" on x's side and it is a 
permutation of the d basis states of the coin space Hd, and the coin C = Cq® In where 7jv acts as 
identity on Hn o,nd Co is a "coin- flip" acting on Hd 

1 

Co = 2\s){s\- Id where |s) = ^V|i). (2) 

^ 1=1 

For a given i S permutes the vertices of the graph, hence 5 is a unitary operation. The permutation 
TT allows us to specify shift operations that act differently on the coin space Hd- Note that the 
coin is symmetric in that it treats all d directions equally, and among all such coins Co is the one 
farthest away from identity. 

Remark: The uniform superposition |$o) = Yli=i 'l2x=i ® 1^) is an eigenvector of U with 
eigenvalue 1 (J7|$o) = l^^^o)); if we start the walk in |$o) it will never change this state. 

To introduce a marked item in the graph we need to have an inhomogeneity in the quantum 
walk by using the coin to "mark" a vertex v, which gives rise to the following: 

Definition 2 [Perturbed Quantum Walk:] The perturbed walk with marked vertex v and "marking 
coin" Ci = —Id is U' = S ■ C, where 

C = Co ® (/ - \v){v\) + Ci ® \v){v\ = C - (Co - Ci) ® \v){v\. (3) 

We will think of U' as the random walk with one (or several) marked coins. This means that instead 
of one coin for all nodes, Co I, we have a different coin Ci on the marked state. Numerical data 
shows that other marked coins exhibit similar properties as Ci = — /, but we will use this coin which 
simplifies the analysis. Then Co — Ci = 2|s)(s|, and U' = U — 2S\s,v){s,v\ = U ■ (idjv — 2|s,u)(s, v|) 
using Co|s) = \s). 

The quantum walk U gives rise to a search algorithm on a graph G in the following way: 



Quantum Walk Search Algorithm 


1. 


Initialise the quantum system in the uni- 




form superposition |$o)- 


2. 


Do T times: Apply the marked walk U' . 


3. 


Measure the position register. 


4. 


Check if the measured vertex is the 




marked item. 



An item on a vertex of the graph could be marked by setting an auxiliary qubit to |1), whereas the 
unmarked items could have this qubit set to |0). Then this auxiliary qubit can control the coin to 
be C for the unmarked items and C' for the marked item. 

We will analyse this algorithm to obtain upper bounds on the query complexity of search by 
random walks. 
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Complete Graph - Grover's Algorithm: As a first example let us illustrate how we can view 
Grover's algorithm |Gro96j as a random walk search algorithm on the complete graph. Each vertex 
has N edges (we will include a self-loop for each vertex). Both vertices and edges are labelled with 
1, . . . , A; the coin space and the vertex Hilbert space are both A-dimensional and we will write 
states as \i) \ j), where the first register is the coin-register. The shift operation S is defined as 

S: \j) \j) ® 

The marked coin in this case is chosen to be Ci = —Co, which gives Ci — Cq = — 2Co and 
C = Cq ® {I — 2\v){v\), where \v) is the marked state. Note that Co = 2|s)(s| — Iat is the 
reflection around the mean operator of Grover's ("standard") algorithms and I — 2\v){v\) =: Ru 
the phase flip of the oracle. Recall that Grover's algorithm is of the form [R^ ■ Co)^|s). The 
initial state for the random walk based algorithm is the uniform superposition |$o) = \s) \s). 
Now C/'|$o) = S ■ C'l^o) = R^\s) Co|s), C • f/'|$o) = (Cq • Ry)\s) {R^ ■ Co)\s) and W^l^o) = 
{R, ■ Co)\s) (Co • So we see that a random walk in this scenario gives exactly Grovers 

algorithm on both the coin space and the vertex space, at the expense of a factor of 2 in the number 
of applications. 



3 Results in 2 dimensions 



We give several upper and lower bounds for the discrete quantum walk on the grid. The A 
memory locations are arranged in a \/A x \/A grid G, labeled by their x and y coordinate as 
\x, y) for x, y G {0, . . . , \/A — 1}. will assume periodic boundary conditions and operate mod \/A. 
The natural coin space is 4-dimensional. We will label the edges emanating from each vertex with 
^) ^5 T) i) indicating the positive and negative x and y directions. 

As it turns out, the choice of the coin transformation (or, equivalently, of the permutation vr in 
Eq. Q) is crucial for the performance of the random walk. We will show that using a "flip-flop" 
shift, gives a search algorithm that succeeds in 0(\/Alog A) time. The "flip-flop" shift Sff changes 
direction after every move, i.e. vr flips | with [ and with <— . Our analysis of the "flip-flop" 
based walk follows the numerical discovery of its performance by Neil Shenvi |She03j . Another 
natural shift is the "moving" shift Sm which does not change direction (i.e. in Eq. ^ tt = id and 
Ivr(i)) = N)). 



Sff-. \^)<^\x,y) — > I ^) (g) |x Sm'- \^)®\x,y) — > \ ^) \x + l,y) 

I ^) ® \x, y) — > \ ^) ®\x - l,y) \^) ® \x, y) — > I ^) |x - 1, y) 

I T) O \x, y) — >\[)®\x,y+l) It)® \x, y) — >\])®\x,y + 1) 

I i) ® \x,y) — >\]) ®\x,y - I) I i) |x,y) — > \ [) ®\x,y - I) 



(4) 



Surprisingly we will show that the "moving" shift gives a walk search algorithm that takes time 
0(A). So even though it seems this walk "moves faster" than the "flip-flop" walk, the resulting 
algorithms performs much worse, no better than classical exhaustive search. It is this surprising 
behavior of Sm which has halted the progress in flnding a good discrete quantum walk search 
algorithm on the grid. 
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Theorem 1 For the quantum walk search algorithm associated to the quantum walk U = Sff ■ C , 
with Sff as in Eq. there is a T = 0{\JN log N), such that after T steps the probability to 

determine the marked state is pT = 0(1/ log A^) . 

Corollary 1 We can get a local search algorithm based on the quantum walk that finds the marked 
state with constant probability in time 0(\/iV log A^). 

Proof of Corollary [Tl The initial state 1$) can be generated with -v/iV local transformations. 
Since we only have an estimate for T up to a constant factor, we need to repeat the random walk 
an appropriate (constant) number of times. For the algorithm we will use amplitude amplification 
|RHMTn2j to achieve a time 0{Vn log N). We will give more details in the proof of Theorem ^ 

■ 

Theorem 2 The quantum walk search algorithm associated with Sm as in Eq. ^ takes at least 
Q{N) steps to determine the marked state with constant probability. 

We also consider a two dimensional coin inspired by Dirac's equation in 2+1 dimensions. Let 
I ff) = |0) and I JJ.) = |1) be the standard basis for one qubit and | <^) = "^|0) + "^|1) and 
I =^>) = ~ "^ll) be the Hadamard basis. If there is no marked coins, one step of the quantum 

walk U with the two-dimensional coin consists of: 

1. Move up /down: 

I t) ^ \x) (S) \y) ^ \ i[) \x)\y-l), 
I -11-) \x) |y) ^ I J|) (8) \x)\y + l). 

I ^) \x) (g) \y) ^ \ (E) \x - l)\y), 
I ^) \x) |y) ^ I ^) |x+ l)\y). 

If there is a marked coin \v), we define the quantum walk as U' = U{I — 2\s,v){s,v\) where U is 
the walk with no marked coin and |s) is the state "^1 fT) + "^1 

Theorem 3 The associated quantum walk search algorithm takes 0{\/N log A^) steps and the prob- 
ability to measure the marked state is 17(1/ log A^). This yields a local search algorithm running in 
time 0{Vn log N). 



2. Move left /right: 



4 Results in 3 and more dimensions 



In more than 2 dimensions the "flip-flop" based quantum walk search algorithms achieves its optimal 
performance of 0(\/iV)- Here G is a grid of N vertices, arranged as x . . . x with periodic 
boundary conditions, as before, and states are labelled as . . . ,Xd). 

Theorem 4 Let G be the d-dimensional grid with N vertices. Then the associated quantum walk 
with one marked coin takes 0{^/N) steps and the probability to measure the marked state is constant. 

Theorem 5 The results of Theorems^ [3 and^ hold also for two marked items. 
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5 Abstract search algorithm 



Before giving the technical details let us give some intuition of the proof. Recall that Grover's 
algorithm in its standard form is a succession of reflections i?^ around the marked state \v) followed 
by a reflection around the mean i?|<j,^ = 2|<I>)(<I>| — Ijsr, where |<&) is the uniform superposition over 
all items. It can be viewed as a rotation in a two dimensional space, spanned by the marked state 
\v) and the initial state. In the basis where |0) = |<I>) and \v) = "^|0) + ^J^^\l), Grover's 
algorithms corresponds to the transformation (with sin0 = 2 ^^~^ ) 



U= "^^'"^ "^^^^'^ . (5) 




The two eigenvectors of U are | it w) = "^(|0) =bi|l)) with eigenvalues e^**^. The initial state is a 
uniform superposition of the two eigenvectors \^) = "^d"^) + I ~ ^))- After T applications of U, 
with T chosen such that Tip = ^ , we have 

u^'m = f/'^^d^) + 1 - ^)) = ^(-^\^) + ^1 - ^)) = ii) 



which has an overlap of y ^^^^ with the marked state | v) . 

In the random walk algorithm the transformation (/— 2|s, v){s, i"!) is a counterpart of and the 
transformation U is an "imperfect" counterpart of We will first show, that with an appropriate 
choice of coin (as in Thms. Q 01 and |1J) the resulting transformation is still approximately in a 
2-dimensional subspace; In this space U' will correspond to a rotation as in Eq. ((SJ- Chosing T 
appropriately will (approximately) give a state with a "large" overlap with the marked state or its 
neighbors. 

In the case of the "bad" coin, as in Theorem [21 we will show that there is a large eigenspace 
of eigenvalue 1 of the perturbed walk, and that the initial state has a large overlap with this 
eigenspace. Hence the state of the system nearly doesn't change by the walk. 

More formally, an abstract search algorithm consists of two unitary transformations Ui and U2 
and two states {ipstart) and {tpgood)- We require the following properties: 

1. Ui = I - 2\tljgood){i'good\ (in other words, Ulipgood) = -\i^good) and, if is orthogonal to 
\Tpgood), then 1-0) = IV')); 

2. U2\ipstart) = \tp start) for some State {ipstart) with real amplitudes and there is no other eigen- 
vector with eigenvalue 1; 

3. U2 is described by a real unitary matrix. 

The abstract search algorithm applies the unitary transformation {U2Ui)'^ to the starting state 
li^start)- We claim that, under certain constraints, its final state {U2Ui)'^\'4! start) has a sufficiently 
large inner product with {ipgood)- 

The next lemmas, which we will prove in Sec. [71 describe the main properties of an abstract 
search algorithm that we use. Let U' = U2U1. Since U2 is a real unitary matrix, its non-itl- 
eigenvalues come in pairs of complex conjugate numbers. Denote them by e^*^^, e^*^™. Let 

dmin = min(0i, . . . , 9m)- 
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Lemma 1 Define the arc A as the set of e'^ for all real 6 satisfying —Omin <d < dmin- Then U' 
has at most two eigenvalues^ in A. 

The two eigenvectors with these eigenvalues will be very important to us. We will show that 
the starting state is close to a linear combination of them. Therefore, we will be able to determine 
the evolution of the starting state by studying these two eigenvectors. 

We start by bounding the two eigenvalues. Let and \^J) be the eigenvectors with eigen- 
values e*^J and e~*^J, respectively. We express IV'good) as a superposition of the eigenvectors of 1/2- 

m 

Lemma 2 It is possible to select |$^) and so that = aJ and a'j is a real number. 

In the next lemmas, we assume that this is the case and denote a'j = aJ simply as aj. 
Lemma 3 The eigenvalues of U' in A are e^*" where 



a = e 



1 



(7) 



Let \wa) and \w-a) be the two eigenvectors with eigenvalues e*" and e respectively. Define 

Wstart) = -^l^a) " ^^-a), \Wstart) = \\u,'^l^J\ l^'start)- We claim that \Wstart) is close to the 

starting state li^start)- This is quantified by the following lemma. 



Lemma 4 A ssume that a < n^min- 

Then, 

{i^start\Wstart) > 1 - 6 (^a^ ^ i^Jl^osO,? 

The last lemma shows that, after repeating U2U1 a certain number of times, the state has 
significant overlap with iV'good)- Say we apply {U2U\)^'^/'^"'^ to the state -^\wa) — -^\w-a)- Then, 
we get the state which is equal to 

-^e^'^/^\Wa) - e~^'^/'^-^\W-a) = ii-^\Wa) + -^Iw-a)) =■ \Wgood) 

plus a state of norm 0{a) (because 7r/4 and [7r/4a]a differ by an amount which is less than a). 
Lemma 5 Assume that a < \6jnin- Let \wgood) = ^l^a) ^~ "T^l^""^' 

\{i^good\wgood) \ = I min I , == I , 1 



^The next lemma implies that there are exactly two eigenvalues in A. 
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Corollary 2 Assume that a < 



2 "mtn ■ 



KV^soodKC/aC/Or^/^^T \w,ood) = e mill , . 1 + 0(a) 



These three lemmas are the basis of our proofs. In each of our positive results, we first find 
a subspace li. such that the search algorithm restricted to this subspace is a special case of an 
abstract search algorithm. Then, we apply Lemma 0] to show that the starting state is close to 
\'WstaTt) and Lemma 131 and corollary [2 to show it evolves to a state having significant overlap with 

6 Proofs of the main results 
6.1 Theorem 1 

Let us determine the eigenspectrum oi U = Sjf ■ (Co In) first. 



Claim 6 [Spectrum of U :] U has eigenvalues X^i with corresponding eigenvectors of the form 
\vH)^\Xk)^\xi) forallk,l = 0,...,VN-l, where \xk) = ^'"'U) with oj = e"^'/^ , 

and Xki and \vki) satisfy the equation 

uj-'' \ 

Ckl\Vkl) = n -I ■ Co\Vkl) = Xkl\Vkl)- (8) 








to ^ 








\ 
































v 













/ 



The four eigenvalues Xki of C^i are 1,-1 and 6^*^*=' where cos O^i = i(cos ^= +cos ■^)- Let \vli), 

l^fci^) ^'^''^ I^h) vectors \vki) for the eigenvalues 1, -1 and e*^^, respectively. Then, is 

orthogonal to \s) for {k,l) ^ (0,0) and is orthogonal to \s) for all {k,l), including (0,0). 

Proof: Apply U to a vector of the form \vki) ^ \Xk) ^ \xi) and note that Sff\ j) \Xk) ^ \xi) = 
uj'^'l i) O \Xk) ® \xi), and Sff\ i) \xk) ® \xi) = '^''1 T) ^ \Xk) ® \xi), and similarly for the y- 
coordinate, which gives Eq. ©. Solving the equation \Cki — XI\ = for A gives the eigenvalues. 
For {k,l) 7^ (0,0) the 1-eigenvector is proportional to (cu'^ [lu'' — 1) , 1 — u>\ io\l — lo^) , lo'' — 1) and 
hence orthogonal to |s) = ^(1, 1, 1, 1). The —1-eigenvector Iv'f^l) is proportional to (w' + 1,uj''{ui'' + 
l),-{io'^ + l),-Lo''{uj^ + 1)) and hence orthogonal to \s) = i(l, 1,1,1). For {k,l) = (0,0), \s) is a 
1-eigenvector of Cqq. ■ 

For {k,l) = (0,0), the eigenvalue 1 occurs 3 times. Thus, there is a 3-dimensional 1-eigenspace. 
Since \s) is orthogonal to \vqq), |s) belongs to this eigenspace. We choose \vqq) = \s) and |u^q) 
orthogonal to \s). 

Let TC'q be the space spanned by the eigenvectors \v^i) (gi \xk) ® \xi)j (kj) ^ (0,0) and |<I>o) = 
boo) ® \Xk) ® \xi)- Notice that all other eigenvectors of U are orthogonal to \s,v), by Claim El 
Therefore, \s,v) is in TC'q. Moreover, applying U' keeps the state in TC'q, as shown by 

Claim 7 We have [/'("Hq) = TC'^. Furthermore, U' has no eigenvector of eigenvalue 1 in TC'q. 
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Proof: For the first part, notice that U' = U{I — 2\s,v){s,v\). Therefore, it suffices to show 
UiTt'o) = TL'q and (I — 2\s,v){s^v\){'H'q) = TL'q. The first equahty is true because TL'q has a basis 
consisting of eigenvectors of U . Each of those eigenvectors gets mapped to a multiple of itself which 
is in TL'q. Therefore, U{H'q) = TL'q. The second equality follows because (/ — 2|s, t;) (s, uDiV') = 

— {s.,v\il))\s,v) . This is a linear combination of \'4j) and \s,v) and, if IV') G TL'q, it is in TLq. 

For the second part assume \ojq) is an eigenvector of eigenvalue 1 of U' in TL'q. Then 

0/ (^ol'^o) = (^o|t/Vo) = {^(i\U{I -2\s,v){s,v\)\ooq) = {^q\loo) -2{^q\s,v){s,v\ooo)- 

This implies {(^q\s,v){s,v\ujq) = and, since {(^q\s,v) = / 0, that {s,v\loq) = and l^o) = 
U'\ijjQ) = U\ojq) which in turn implies that |a;o) is an eigenvector of eigenvalue 1 of U . Since |a;o) has 
zero overlap with \s,v) and precisely the 1-eigenvectors of U orthogonal to |$o) have zero overlap 
with it follows that (<I>o|u;o) = which contradicts that |u;o) £ "^o- * 

The above shows that the random walk algorithm starting in |<I>o) is restricted to a subspace TL'q 
of the Hilbert space. Since is the only 1-eigenvector of U in TL'q, we have an instance of the 
abstract search algorithm on the space TL'q, with Ui = I — 2\s,v){s,v\, U2 = U, \Tpgood) = \s,v) and 

\lpstart) = l^o)- 

As described in Section [3 we study the 2-dimensional subspace spanned by \wa) and \w^a)- 
First, we bound a using Lemma |21 We need to expand {tpgood) = \s,v) in the basis of eigenvectors 
of U. Define = \v^i) <8> \xk) ^ \xi)- Let 1^"^;) be the vector obtained by replacing every 

amplitude in by its conjugate. Then, = \vZk -i) "25 \x-k) "X" |x-/)- (This follows from 

two observations. First, replacing every amplitude by its conjugate in \xk) gives \x-k)- Therefore, 
\^ki) = \v)^\X-k)(^\X-i)- Second, since f/l^^;) = 6*^*='!$^;), we have = 6"*^*=' |#^;), implying 

that \v) = _/)•) From Lemma |3 we have 

\s,v) = ao\<^o) + and^fci) + l^fei)) 

(fc,0^(o,o) 

where |$^;) and |$^;) appear with the same real coefficient aki = {s,v\^^i) = {s,v\^^i). 
Claim 8 ^ 



/2N 

Proof: We have \{v\xk) ® \xi)\ = (since each of the locations has an equal emplitude 
ii^ \Xk) ^ \xi))- It remains to show that = For that, we first notice that |s) is a 

superposition of \v^i) (since and {vj^i^) are orthogonal to \s)). By direct calculation {s\Cki\s) = 
i(cos ^ + cos We have 

This is possible only if = Ks|v^;)| = 73- ■ 

Therefore, 

\s,v) 



1 > 1 



{k,l)f^{0,0) 

By Lemma El a = 0( , ^ ^ ). The following claim implies that a = Q{-^=^^=^). 



kl 1-COS9 
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Claim 9 E(.,/)^(o,o) = @{NlogN). 

Proof of Claim IHJ Recall from Claim |H1 that the eigenvalues corresponding to \v^i) \xk) \xi) 
are e^*^*-' = cos 9ki ± i sin O^i where cos 9ki = |(cos ^= + cos ■^)- For x G [0, 2tt], we have 

l-^<cos2;<l ^. (10) 

2 TT^ 



Therefore, 



±{k' + i^)<cose,,<i--L-^k^ + i^), 

+ /2) < 1 _ cos < + l') (H) 



This means that it suffices to show 



^Efc2^ = 0(^l°g^)' (12) 



where the summation is over all k, I G {0, . . . ,\'N — 1} such that at least one of k, I is non-zero. 
This follows because E^i. ^ pq:^ = 0(logA^). A simple way to see this is to sum points that lie 
on m-rectangles with the four corners (ibm, ibm). The term pqip- for {k,l) on an m-rectangle is 
bounded as < -^^^ < and there are 8m such points on each m-rectangle. Hence 



N-l ^ ^ VN-1 ^ 

m=l k,l m=l 

The claim now follows from X^m^]"^ m~\ ~^ '^(1))- ■ 

Next, we use Lemma|l]to bound the overlap between |<I>o) and \wstart) = '^l'"^") ~ 

Claim 10 E(.,o^(o,o) (railF = ®(^')- 
Proof of Claim ITOl : Using the proof of Claim IHl 

7r4iV2 1 4A^2 



< -. < 



Therefore, it suffices to bound N'^'}2{ki)^{0Q) [k'^+i'iyi ■ Again, we sum points {k,l) over rectangles 
with corners (ibm, ibm). Each rectangle has 8m points, each of which contributes a term of order 
-K to the sum. Since 8m^ = is bounded by a constant, the lemma follows. ■ 

This means that the overlap between the starting state and \wstart) is 1 — Q{a^N'^) = 1 — 
0(i^p-^)- Equivalently, |^>o) = \w start) + rem) , with |||^>r-em)|| = ©(iHilv)- Aft^^ repetitions, 
the state becomes \wgood) + Wrem)^ with H'l'^emll = ®(lHglv) + 0{a) = ©(y^^J-^v)- Finally, we bound 
{wgood\s-,v) , using Lemmaini Since all au are equal to -^7= and cotx < ^, we have 

1 1 

> „ . • (14) 
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From the proof of Claim |^ we know that is bounded from below and above by j 



const 



-cos ffej 



Therefore, Claim El implies Y.kiW = Q{N\ogN). Thus, the expression of Eq. ()14l) is of order 

_^ ' kl 

To conclude the proof of the theorem, the overlap of the state of the algorithm after [^] steps 
and |s, v) is 

\{Wgood\s,v) + > \{Wgood\s,v)\ - \{^'^^^\S,V)\. 

The first term is of order Vl[-j^=) and the second term is of lower order (^( i^g jy ))- Therefore, 
the result is of order Vl( ., .. ). 

^ Vlog N ' 

Hence a measurement gives the marked location \v) with probability p > 1^5^77 • This completes 
the proof of Theorem ^ 



Proof of Corollary 1: First, note that it is possible to generate the initial state |<&o) with 2\/ N 
local transformations. We start with the state concentrated in one point (say |0, 0,0)) and first 
"spread" the amplitude along the x-axis in \fN steps. In the first step we rotate the coin register to 
-3^|0) + ^^^^^ 1-*^)' followed by a |l)-controlled shift in the x-direction, followed by a rotation of 
the coin register back to |0) in the vertex (0, 0). Similarly we repeat this procedure to move ^J'^^^ 

of amplitude from (1, 0) to (2, 0) and so on. After y/N steps we have a uniform superposition over all 
vertices with y-coordinate 0. We repeat this process for the y-direction, which gives us the uniform 
superposition after another yfN steps. Note that this procedure also allows us to implement the 
reflection around the mean, = I — 2|<I>o)(*l?o| in 4\/^ steps: we simply run the procedure in 

reverse (which maps to |0,0,0)), then invert the state |0,0,0) (which can be done locally in 
the vertex (0,0)), and run the procedure forward again. 

Note that we have determined the run-time T only up to a constant (using Eqs. Ullj) and 
(El, (ESI) we can bound Tmin = ^^^^^ <T< ^V^j|s^ = Tmax)- To get e-close to the state 
U''^\^o) we use a standard trick and run the walk for times Tmin, (1 + £)Tmin, (1 + £)'^Tmin, ■ ■ ■ 
until we reach Tmax- One of these times is within a factor of (lie) of T and hence our state 
and final measurement probability will be e-close to the state at time T. We can chose e to 
be some small constant. The total time including all repetitions (bounded by j^Tmin) is still 
0(T) = 0{^N\ogN). 

Finally, to amplify the success probability we will use amplitude amplification |BHMT()2] , which 
is a succession of steps consisting of reflection around the mean |<I>o) and a run of the algorithm. 
The intermediate reflection around the mean can be implemented in 4\/]V steps, the random walk 
takes 0{\JN log A^) steps, and we need 0(\/Iog7V) rounds of amplification to obtain a constant 
probability of success, which gives a total running time of 0{-\fN \ogN). ■ 



6.2 Theorem 2 

Proof: The key difference between this walk using Sm and the walk from Theorem ^ using Sff 
is that the initial state |<I>o) now has very large overlap with the eigenspace of eigenvalue 1 of C/ 
and U' . This means that the walk (nearly) does not move at all and the state at any time T has 
overlap with |$o) close to 1. The difference becames apparent in the eigenspectrum of U: 
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Claim 6' [Spectrum of U:]U has eigenvalues Xm with corresponding eigenvectors of the form \vki)'Si 
\Xk) \xi) forallk,l = 0,...,VN- 1, where \xk) = ^ E}2"' ^'"'Ij) ^'^^ ^ = e^^'/^ , and \m 
and \vki) satisfy the equation 



Ckl\Vkl) 
















\ 

























J 







v 













) 



Co\vki) = \ki\vki) (15) 



The four eigenvalues Xm of Cm are 1,-1 and 6^*^*=', where cos O^i = — |(cos + cos ■^)- For 
the eigenvector corresponding to eigenvalue 1, we have 



\{vli\s)\> 



1 I 2-Kk I 2ttI I 27r(A;+0 
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Proof: The first part is by straightforward calculation as before (Claim EJ. For the second part, 
the eigenvector corresponding to eigenvalue 1 is If^^) = |j~T^ with 

14;) = {w\l + w^),l + w\w^{l + w^),l + w'') 

We have < 4 because each of the 4 components of |u^;) is at most 2 in absolute value. It 

remains to bound (n^Js). We have 

{vli\s) > = -{w'^il + w^) + 1 + + w^il + w'') + 1 + w'') = -(1 + w'')il + w^). 

4 8 4 

The real part of this expression is 5(1 + cos ^= + cos + cos ^^^^'^ ) ■ This implies the claim. ■ 
Let Hi be the 1-eigenspace of U, spanned by the |<I>^;) = ^ \xk) ^ \xi) iov k,l = 0, . . . , \fN — 1, 

s,v) = "^aki\^li) + \s') 

k,l 

where \s') has no overlap with the 1-eigenspace TCi. We claim 



with |$qq) = |<I>o). Write 



Claim 11 U' has a 1-eigenvector |$-^) such that \{^o\^ )\ =1 • 

Proof: Let Pki = 1 Let |<I>) = / be the projection of \s,v) on Tii. Since 

(<l>o|*^*) = Poo, we can write 

l^-o) = /?oo|^) + Vl- l/?oop|^^> 

where \^~^) is a vector perpendicular to |$). Since |$o) and |^) are both in the subspace Tii, 
l^"*") is also in Tii. We claim that l^"*") is a 1-eigenvector of U'. The state l^"*") is orthogonal to 
\s,v) because l'^-'-) belongs to Hi and is orthogonal to |<I>) which is the projection of \s,v) to that 
subspace. Therefore, 1$"*") is a 1-eigenvector of I — \s,v){s,v\. I^*"*") is also a 1-eigenvector of U 
because it belongs to Hi. This means that it is a 1-eigenvector of U' = U{I — \s, v){s, v\) as well. ■ 
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To complete the proof, we need to bound |aoo|i l^oili • • - i |« v^-i Vjv-iI' have aoo = We 
wiU show that there are Q{N) other a^i of order Q{1/\/N). This would imply that the overlap of 
|$o) with a 1-eigenvector of U' is 

- 1 - f^(^). 



Claim ini gives the desired a bound on the Uki- 

\(s\v^)\ 1 + cos^ + cos^ + cos^^^ 
K^,^bi/) ^ IXfc) ^ 1X^)1 = \{v\Xk ^ Xi)\ X Ks|^;i,)| = > v^F 



4ViV 

e range for ^= and ^= is [— f,f]- Therefore, for half of all k (resp. half of all /) we have 
< f (resp. < f )• For the ^ pairs (A;,/) that satisfy both of those conditions we have 



2Trk 2-kI 2TT(k + l) 11 ^ 

1 + cos —= + cos —= + cos — ^ ^ ^ >l + ^ + ^>l + V2. 
\/iV //V //V \/2 a/2 

Thus, for at least ^ pairs (A;, /) 

1 + V2 1 



4VN 2\/]v' 



6.3 Theorem 3 

Proof: The proof for this random walk algorithm with a 2-dimensional coin proceeds in close 
analogy to the proof of Theorem^ and we will emphasize and prove the points that differ. 



Claim &^ [Spectrum of U:]U has eigenvalues A^^ with corresponding eigenvectors of the form \Vf^i)(^ 



\Xk) ® \xi) for allk,l = 0,...,^/N-l, where \xk) = ^ E]2"' with uo = e^^/v^^ and \f 



and \v^i) satisfy the equation 

I Jcosk iuj-'-sink V \ , , \ /i^^ 
Ckl\vkl) = . ; . , -I , \vkl) = >^kl\vkl)- (16) 
\ iLO sm k CO cos k I 

The two eigenvalues A^^ of Cki are 6^*^*=' where cos 9ki = ^(cos ^^^^^'^ + cos ^"^^^^ ). 

As a corollary, we have that there are exactly two eigenvectors with eigenvalue 1, both of them 
of the form I^^q) |xo) ^ \xo)- Since the coin space is 2-dimensional, the two vectors |uoo) span it 
and, therefore, \v) \xo) \xo) is an eigenvector for any \v). In particular, we can take \vqq) = \s) 
and \vqq) = js-*-) where Is-*-) _L \s). Similarly to Theorem 1, let Hq be the space \vqq) (8) |xo) ^ \xo) 
and \v^i) (d) \xk) \xi)- Similarly to ClaimEl Ti.Q is mapped to itself by U' . 

Let 1$^;) = (g) \xk) ® \xi) and |<I>^;) = \vZi,_i) iX> \x~k) 'i^ \x~i)- We can express the state 
|s, v) as 

\s,v) = aol^o) + Yl + I'J'h))- (17) 



14 



where the equality of the coefficients of and l^i^i) follows from the proof of Lemma [21 (and 
and 1^^;) being complex conjugates). 
We now have to bound the sum of Eq. (jT)). We claim that replacing all by ^ does not 

change the value of the sum. To see that, we first notice that O^i = O^^^^i. Therefore, replacing 
\a\i\ and |a^^_;| by l^^'l +l°~fc,-il (^Qgg j^q^ change the sum. Futhermore, la^p + \a-k-i\'^ = jq- 
(We have |om| = \a^k-i\ = {^\_i\s,v) = (^'Ifc_Js,u). The vectors |$^;) and I^I^.^) 

are the same as \v^i) (8) \xk) ® \xi)- Therefore, \aki\^ + \a-k-i\'^ equals the squared projection of 
|s, v) to \v^j) ® \xk) ® \xi) which is equal to |(f |xfc) ® lx«)P = jj-) Also, we still have oq = and 
|aoP = jj- Therefore, Eq. simplifies to 




just as in the proof of Theorem ^ We get a = Q{\/N log N) as in Lemma 121 The other two parts 
of Theorem ^ also follow in a similar way. ■ 

Similarly to Corollary 1 we can get a random walk based search algorith that determines the marked 
state with constant probability in time 0(^/1^ log N). 



6.4 Theorem 4 



Proof: Let us call the 2d directions on the d-dimensional grid i±, where i indicates the dimension 
and lb the direction of the walk. Then the "flip-flop" shift operation takes the form 



SffH 



\Xl . . . Xi . . . Xfi) 



\XI . . . Xi ± 1 . ..Xd) 



Let us recapitulate what the key elements in the proof of Theorem ^ are and how they generalize 
to the d-dimensional case. 



Claim &' [Spectrum of U:] U has eigenvalues Xkik2...kd with corresponding eigenvectors of the form 
\vkik2...ka) ® IXfci) ® IXfca) ® • • • ® \Xka) (h = 0, . . . , </N - I), where \xi) = T,k=o~ ^1^) '^^^^ 
u! = e^'^Vv^^ Qjj^ \kik2...ka and \vkj^k2...ka) satisfy the equation 





( 













\ 























Ckik2...ki\vkik2...ki) — 

















• CQ\vk^k2...ki) = ^kik2 























\ 










/ 





The 2d eigenvalues \kik2...kd of Ck-ik2...kd '^^^ ^"-^ ~^ '^^^^ multiplicity d — 1 each, and e^*^''i'=2- -fed 
where cos6kj^k2...kd = ^Yli=i'^'^^^^- I'^fcifc2 kj corresponding to eigenvalue 1 are orthogonal 
to \s) for {kik2 ■ ■ ■ kd) 7^ (0, 0, . . . , 0). All \v[^^j^^ f^J corresponding to eigenvalue —1 are orthogonal 
to \s) for all {kik2 ■ ■ ■ k^)- 
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Proof: Eq. H18() is obtained in the same way as in the proof of Claim El Ckik2...ka consists 
of 2 parts, the 2-block-diagonal matrix (call it Dkj^k2...kd) ^^'^ ^o- Each block in Dk^k2...kd ^^s 
eigenvalues ±1 and eigenvectors , itu;^), so the matrix Dk^k2...ka itself has eigenspaces of 

1 and —1 of dimension d each. In each of these two eigenspaces we can find d — 1 orthogonal 
vectors orthogonal to \s). For those vectors Co = 2|s)(s| — / just flips their sign, so the — 1 
eigenvectors of -DA:ifc2...fcd become +1 eigenvectors of Ckik2...kd ^^^^ versa. Call the remaining 
two eigenvectors of eigenvalue ±1 \e±) and expand \s) = Q+|e+) +a_|e_). Then {s\Ckik2...kd\^) — 
{s\Dk^k2...kd\s) = ^Xlf=iCos27r/cj/v^ = cos 6*^1^2... fcd implies |a+p - |a_p = cos6kik2...kd- Let 
be an eigenvector of Ckik2...kd "^^^ orthogonal to \s) and expand = /3+|e+) + /3„|e^). Then we 
obtain (using |e+) = a^|s) + a-\s-^) and |e_) = q*_\s) — a+ls"*") where Is"*") is orthogonal to \s) in 
the space spanned by \e±)) for the eigenvalue {uj\Ckik2...kdW) — (|a+p — |a-P) + 4iIm/3i!^/9_a+Q;!L = 

COs6'fciA;2...A:d =t ^Sin6'fcifc2...fcd- ■ 

For {kik2 ■ ■ ■ kd) = (0, 0, . . . , 0) there are d + 1 1-eigenvectors, we set |fo...o) = N). Then, the 
other 1-eigenvectors are orthogonal to \s). 

Similarly to Theorem^ we restrict to the subspace TC'q spanned by \v^^j^^ (X" Ixfci) "SD \ Xk2) 
■ ■ - ^ \Xkd) and |foo...o) ® Ixo) <X Ixo) ^ ■ ■ - ^ \xo)- As in TheoremQ] we have U'{Hq) = Hq. Further 

\S,V) = ^|$o) + J2 (\Kk2...kd) + \^k,k2...kd))^ 

^ ^ {ki,k2,-,kdm0,0,-,0) 

where 1$^^^^.,.^^) = \vl\,,„k,) ^ IXk^) ® IXfea) ® • • • ® \Xkd) and 1$^^^^,..^^) = \v^X2...kd) ® l^-^i) ® 
\X-k2) (8) • • • (X" \x-kd)- We use a modified ClaimlHJ 

Claim 9' E(fc„fc2,...fcd)^{o,o,...,o) i^cosel^,^...,^ = 

Proof: The proof follows along the lines of the proof of Claim IHl By Claim IHI", 



l-cos0fc,fc2...fcd E?=i(l 



2TTkj < 

COS -^TJd, 



1 



Similarly to LemmaEl this is bounded from above and below by a constant times N'^^'^ ,-2.,-2. .i,-2 - 
Thus, we have to estimate 

N'/' E ,2 ,2 ^ 72 (19) 

where the summation is over all ki S {0, . . . , ~ 1} such that at least one of the ki is non- 
zero. We divide tuples the (ki, . . . ,kfi) into N^/'^ "slices", with the m**^ "slice" containing those 
tuples where max(A;i, A;2) • • • 1 ^d) = w-- The m*^ slice contains 0{m'^~^) tuples. Therefore, the 
sum y^. . , » , , — -T7 over the m^^ slice is of order m^~^-^ = m'^~^. Since there are N^/'^ 

^ki,...,kd kf+k^+...+k^ 

slices and, for each of them, the sum is m'^~^ < N^'^~^^^'^, the sum k^+k^+ +k^ 

(/ci, . . . , fed) / (0, . . . , 0) is of order at most iVi/<iAr{'^-3)/d ^ j^(d-2)/d^ -pj^jg implies that Eq. ^ 
is of order at most N . It is also of order at least since each individual term N'^^'^ i,2,i.-i\ .u-i is 
at least a constant. ■ 
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Therefore, applying Lemma 01 gives a sharper bound a = Q{-^). Notice that we still fulfill the 
requirement a < ^Omin needed for Lemmas H] and [3 The reason for that is that all Oi are at least 

Similarly to Claim EJ, we can show 

(A;i,A:2,.--A;d)7^(0,0,...,0) i ^ «/ 



and 



(fci,fc2,.-fcd)7^(0,0,-,0) kik2...kd 



By combining these two equalities with Lemmas ^ and El we get that the overlap between the 
starting state |$o) and \wstart) = ■^\wa) — ^\w-a) is 1 — 0(-^) and the overlap between \wgood) = 
-^\wa) + ■^\w-a) and \s,v) is This implies that the search algorithm's final state has a 

constant overlap with \s,v). ■ 

6.5 Theorem 5 

We first discuss generalizing the positive results (Theorems ^ Bl and (HI) to the case with two 
marked items. The main issue is to state them as instances of the abstract search. Assume 
there are k marked locations vi,...,Vk- Then, one step of the search algorithm is U' = {I — 
'2J2i=i \s,Vi){s,Vi\)U. Currently, we are not able to analyze cases of the abstract search where U2 
flips the sign on more than a 1-dimensional subspace. 

For the k = 2 case, we can avoid this problem, via a reduction to A; = 1. Define \s') = 
-^|s,ui) + -^\s,V2)- We claim that applying {U')^ to the starting state |$o) gives the same final 
state as applying (C/")^ where U" = {I - 2\s'){s'\)U. 

To show that, let T be a symmetry of the grid such that T{vi) = V2 and T[v2) = vi. (For the 
2-dimensional grid, if vi = (xi,yi) and V2 = {x2,y2), then T(x,y) = (xi + X2 — x,yi + y2 — y)-) We 
identify T with the unitary mapping |c, f) to |c, T(t>)). 

Claim 12 For any t > 0, r([/')*|$o) = (U'Yl^o)- 

Proof: By induction. For the base case, we have to show T|<I>o) = I'^'o)- This follows since |$o) is 
a uniform superposition of the states \s,v) and T just permutes locations v. 

For the inductive case, notice that T commutes with both / — 2|s, vi){s, vi\ — 2\s, V2){s, V2\ and 
U. Therefore, TU' = U'T. If the inductive assumption T(U'Y\^o) = (f/')*|^>o) is true, then we 
also have 

T([/')*+^|^'o) = U'T{U'Y\<^o) = {U'Y+^\<^o), 
completing the induction step. H 

Let \s") = -^\s,vi) - -^\s,V2). Then, (C/')*|^o) is orthogonal to \s") because T\s") = -\s"). 
We have \s, vi){s, vi\ + \s, V2){s, V2\ = \s'){s'\ + \s"){s"\. Together with {U'Y\^o) _L \s"), this means 

(I -2\s,vi){s,v,\ -2\s,V2){s,V2\){U'Y\^o) = {I-2\s'){s'\){U'Y\^o)- 
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Thus, U' can be replaced by U" at every step. 

The rest of the proofs is now similar to the case of 1 marked location, except that \s,v) is 
replaced by \s') everywhere. Theorem |21 also follows similarly to the 1 marked item case, with \s') 
instead of \s,v). 

7 Proofs of the technical lemmas 

In this section, we prove Lemmas ^ 13 Hi and [3 We will repeatedly use the following result which 
can be found in many linear algebra textbooks. 

Fact 1 The eigenvectors of a real unitary matrix either have eigenvalue ±1 or else they appear in 
conjugated pairs with eigenvalues e^*"^ and eigenvector | it a;) = -^{\R) i i\I)), where \R) and \I) 
are real normalised vectors and {R\I) = 0. 

Proof of Lemma [TJ Let g = I — i?e(e*^'"™). Then, for any e*^ G A, Re{e^^) > I — g and, for 
every eigenvector |a;) with an eigenvalue in A, Re{uj\U'\Lj) > 1 — g. 

If there were more than two eigenvectors of U' in TC'q with eigenvalues on the arc A, we 
could construct a linear combination \a) of them such that \a) _L |<I>o)) N, w). Since \a) is a lin- 
ear combination of vectors \uj) with Re{uj\U'\uj) > 1 — g we have Re{a\U'\a) > 1 — g. But then 
Re{a\U'\a) = Re{a\U{I — 2\s,v){s,v\)\a) = Re{a\U\a). Since \a) is orthogonal to |<I>o) and all other 
1-eigenvectors of U (|a) G Wq), \a) is a linear combination of eigenvectors of U at least g away from 
1 and hence Re{a\U\a) < 1 — g, which gives a contradiction. ■ 

Proof of Lemma 121 Let |<I>) be the vector obtained by replacing every amplitude in \^^) by its 
complex conjugate. Since U2 is a real unitary matrix, f72|^j^) = e*^j|<&^) implies f/2|^) = e~*^j|<^). 
Therefore, we can assume that |<I>~) is a complex conjugate of The coefficients and aj 

are equal to the inner products start) and {^J\ipstart)- Since \ipstart) is a real vector, these 

two inner products are complex conjugates and = (oj)*- By multiplying and \^J) with 
appropriate constants, we can achieve a'j' = aJ . ■ 

Proof of Lemma 151 First, we express \ipgood) in the basis consisting of eigenvectors of U2- 

m 

\^good) = aol^fo) + E^J-d*^/) + l*^/))- (20) 
i=i 

where |$o) = start)- We define for real a 

10 =aocot-|<i>o) + 2j«j (cot^-^|<I>+) + cot^-^|cI>-)j . (21) 

i 

Similarly to Claim 2 in |Amb03j . we have 

Lemma 13 If \w'^) is orthogonal to {ipgood), then \uja) = \ipgood) + A'^'a) eigenvector of U' 

with eigenvalue e*" and \iO-a) = iV'good) + i\w'-a) is an eigenvector of U' with eigenvalue e~*". 
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Proof: The proof is similar to |Amb03j . but we include it for completeness. 
Apply U' to \iVa) and expand in the eigenbasis of U: 



Y^^a^ (e*e.(-l + icot^^^)|$+)+e-^^^(-l + icot^)|^>: 



U'\uJa) = U{I -2\s,v){s,v\){\il:good) + i\w'^)) = U{-\iPgood)+i\w'a)) = aQ{-l + icoi%)\^o) + 

In this equation, every coefficient is equal to the corresponding coefficient in e^^dV'good) + 
Namely, for the coefficient of |<l*o)) we have 

-l + icot-) = — ^ = e*-— ^ = e*-(l + icot- 
z / sm ^ sm ^ V 2 

For the coefficient of \^^)-, we have 

■a I — c-,-+a\ a e^^^~~^^^^ ■ e^^\^~k~^^ ( —B^ + a 



e*''^ -1 + icot ' = ^— = e*" = e*" 1 + icot 

\ ^ J sm — ^ — sm — ^ — 

and, similarly, the conditions for the coefficients of are satisfied. H 
By Eqs. and dH)), {s,v\w'^ = is equivalent to 

flo cot - + ^ a| (cot + cot ) = 0. (22) 

Let 

Smin be the smallest of ^i, . . ., Om- Then, this equation has exactly one solution in [0,^mm] 
and one solution in [ — ^mm^O]* The reason for that is that the cot function is decreasing (except 
for X = kir, where it goes to — oo for x < kir and +oo for x > kir). Therefore, the whole right hand 
side is decreasing, except if one of ^, ° ^ ^ i ~~2^ becomes a multiple of tt. This happens for a = 
and a = zizOmin- Since Omin is the smallest of 6j, [—Omin,0] and [0, contain no values of a for 
which one of the cot becomes infinity. On the interval [0, Omin] the left-hand side of (|22|) goes to 
+00 if q; ^ 0, —00 if a — > ^oi and is for exactly one value of a between and ^oi- This means 
that the two eigenvectors of U' in the arc A are of the form \s,v) + i\u;'^) with \u)'^) as in Eq. (|^T|) . 

Next, let us determine this a. Since 

cos X cos y cos x sin y + cos y sin x sm(x + y) sm(x + y) 

cot x + cot y = h = : : = : = 2 — ; — -, 

smx smy smxsmy smxsmy cos(a; — 7/j — cos(x + yj 

Eq. ll^ is equivalent to a?, cot § + y",- 2a^— — — = which, in turn, is equivalent to 

^^^^^^ ^ ^ * — J C'Oo w j COS CXr 

cotf 1 

sm a ^-^ ■' cos a — cos 

3 ■' 

For Q = 0(1), sina = (1 — o(1))q and coto = (1 + o(l))^. Therefore, we have, with cosq < 1 

(1 + «(i))4 = E ^ ^ E «r T-^- (23) 

^-^ ■' cos a — cos Oj ^-^ ■' 1 — cos Oj 
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This implies 

Q < — (24) 



^3 l-cose, 



It remains to lower bound a. 

Assume that a < ^9min- (Otherwise, the lower bound of the lemma is true.) Then, we have 



cos a — cos 9j > cos — cos 9j > -(1 — cos 9j). (25) 



The first inequality follows from a < < and cos being decreasing on [0,7r]. The second 
inequality is equivalent to cos > ^(1 + cos 9j) which follows from 1 + cos 9j = 2 cos 9j < 2 cos 9j. 
Eq. ((231) and ((23) together imply 

which implies the lower bound on a. ■ 
Proof of Lemma (IJ We will show that the starting state is close to the state \ w start) = tt-^ — ill^s 



start! — |UI7 JT I '"start/' 

II startW 

\w'start) = T^l^c) - By Eq. ((211, we have 



start! 



\/2ao cot '^i\tpstart) + X] (cot "^^^ - cot ^ i - 



We have {^start\w start) = ^^^'^^^^^^^^ = ^J'"' f ■ Therefore, we need to bound \\w',,,,,\\. We 

II startll \\ StartW 

have 

2 



start I 



2al cot^ - + 4 a2 (cot - cot ) 



Since cotx = (1 + o(l))^, the first term is 2(1 + o(l))aQ^ = ©(ag/a^). Similarly to the previous 
lemma, we have 

a + 9j 9j — a sina (1 + o(1))q 

cot — — ^ — cot — — 



cos a — cos Uj cos a — cos 
Since a < ^9min, we have cos a — cos^j > ^(1 — cos6'j) (similarly to Eq. ([25]) ). Therefore, 



2 



a^cot^f^"' G(a2/a2) ^ ^ al \i-cos9j 



This means that 



l|w|l I Og (1 — COSt^j)^ 
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Proof of Lemma [SJ Let \wgood) = "^1^") ~ We consider the unnormalized state 



\<ood) 



good! 



\ - !<) - \w'_^). Obviously, \wgood) = hF^- We have 



\w'good) = 2ao|V'siari) + 

a, (^(2 + icot^-^ + zcot _^)|^+) + (2 + ^cot^-^ + ^cot ^)|^7)J . 



Also {iPgood\wgood) = ^^vj^'^ - Futhermore, {%ood\w'J = 2^0 + Eti = 2|||V'aood>f = 2. 

II good'' ^ J J 

(The imaginary terms cancel out because cot ^"^^^^ = — cot •) It remains to bound H^t'goodll- 

We have 

2 



Since a < ^9min, this sum is at most 



''aoodf = 2ag + ^ + ^ 20^ ( cot + ^ cot ) 

2S (^cot^-^ + cot^-^j • 



2 + 5^ 2a2(2 cot < 2 + e cot 



2 

4 



Therefore, = 0(max(y a?- cot^ ^, 1)) and (ipgoodlw) = 6 j min j , 1 



8 General graphs 

The approach and methods we have presented are amenable to analyze quantum walk algorithms 
on other graphs G. All we need is the eigenspectrum of the unperturbed walk U, an appropriate 
subspace H'q containing no 1-eigenvectors of U (which is equivalent to proving that all but one 1- 
eigenvector of U is orthogonal to \s)), and the sums in Lemmas El 13 and El involving the eigenvalues 
of U (which give the angle a and the overlaps). Then we can apply Lemmas l3l5l to get the desired 
result. 



Hypercube: For instance we can derive the performance of the random walk search algorithm 
on the hypercube, given in jSKWDSj without having to guess the form of the eigenvalues | it ujq). 
|SKWn3] showed that the random walk search algorithm after time T = gives a probability 

of ~ ^ to measure the marked state. For the d-dimensional hypercube (with N = 2'^ vertices), the 
transformation U has d2'^ eigenvectors. An argument similar to Claim [7| shows that the quantum 
walk stays in the 2'^-dimensional subspace spanned by 2'^ eigenvectors with eigenvalues 6*^*= |MR02j 
with 

2k 

cos 9k = 1 - -r 
a 
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for k = . . . d, each with degeneracy (^). Among those, |^>o) is the only eigenvector with eigenvalue 
1, thus we have an instance of the abstract search. We can now apply Lemmas I3I5I For Lemma El 
we need to compute the sum of inverse gaps of all eigenvalues of U in Lemma 1^1 which is now of 
the form 

k=i ^ ^ " k=i ^ ^ 

This gives a time T = 0(\/iV) to rotate the state. Evaluating the quantities in Lemmas ^ and [5] 
shows that measuring the final state gives a marked location with constant probability. 

Remarks on general graphs We have seen that crucial to "good" performance of these algo- 
rithms are essentially two ingredients: 

1. Coin property: The relevant Hilbert space TC'q of the perturbed walk U' does not contain 1 
eigenvectors (i.e. all 1-eigenvectors of U except the starting state \^o) are orthogonal to the 
marked state \s,v)). 

2. Graph property: The gap g between the 1-eigenvalue and the real part of the next closest 
eigenvalue of U is sufficiently large (determines the overlap of the initial state with the two 
relevant eigenvectors). Furthermore the sum of inverse gaps of the eigenvalues of U, i.e. 
of terms (1 — cos6)~^ where 6 is an eigenvalue of U, is sufficiently small, such that the 
angle between the two eigenvectors of U' with eigenvalue closest to one is sufficiently large 
(determines the speed of the algorithm). 

We call the first item a "Coin" property because, as we have seen, it is the choice of coin that 
determines this behavior. We call the second property a "Graph" property because the gap and 
the closeness of the perturbed eigenvalues to 1 depend on the topology of the graph. 

To be more precise let us carry our argument through for Cayley- graphs of Abelian groups. 
The eigenvalues of the unperturbed walk U are "split" by the coin to be "around" the eigenvalues 
of the normalized adjacency matrix of the graph. More precisely, note that the adjacency matrix 
of a Cayley graph can be written as a sum of commuting shift operations over all d directions 
A = 2 Z]f=i ^i- The eigenvalues of A are just sums of eigenvalues of shifts (which are the Fourier 
coefficients). For instance in the case of the grid the eigenvalues of A are of the form ^(cos + 
cos ■^) = \{oj^ + oj~^ + cj' + w^') for A;, Z = . . . ^/N — 1. From Eq. ^ we see that in general 
the coin "changes" the eigenvalues to be some linear combination of oj^^'K This is what happens 
in general, and the resulting eigenspectrum will have gaps of the same order of magnitude as the 
spectrum of the matrix A (which defines the simple random walk on the graph). Thus, it is likely 
that the sum iJ^^ose eigenvalue of U) which determine the speed of the algorithm can be 
estimated from graph properties alone (given a successful choice of coin). 

Related work: Analyzing quantum walks on general graphs has been recently considered by 
Szegedy |Sze04j who has shown a following general result. If e fraction of all vertices of a graph is 
marked and all eigenvalues of a graph differ from 1 by at least 5, then 0(--^) steps of quantum 
walk suffice. This contributes to the same goal as our paper: developing general tools for analyzing 
quantum walks and using them in quantum algorithms. 
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The power of two methods (ours and |Sze04p seems to be incomparable. The strength of our 
method is that it is able to exploit a finer structure of the graph. For example, consider the 2- 
dimensional grid with a unique marked item. Applying the theorem of |Sze04j gives a running 
time of 0{N^/^)^ compared to 0{-\/N\ogN) for ours. (One marked item is a 1/N fraction of all 
items and the eigenvalue gap 5 for the grid is I/^/N.) The reason for this difference is that most 
eigenvalues are far from 1. The approach of |Sze04j uses worst-case (minimum) difference between 
1 and eigenvalues of the graph, which is small {^j^)- Our approach uses a quantity (^g i_cosg ) 
that involves all eigenvalues, capturing the fact that most eigenvalues are not close to 1. 

The strength of Szegedy's |Sze04j analysis is that it allows to handle multiple marked locations 
with no extra effort, which we have not been able to achieve with our approach. It might be 
interesting to combine the methods so that both advantages can be achieved at the same time. 

9 Conclusion 

We have shown that the discrete quantum walk can search the 2D grid in time 0{^/N \ogN) and 
higher dimensional grids in time 0{^fN). This improves over previous search algorithm and shows 
an interesting difference between discrete and continous time quantum walks. More generally, we 
have opened the route to a general analysis of random walks on graphs by providing the necessary 
toolbox. 

The main open problems are applying this toolbox to other problems and learning to analyze 
quantum walks if there are multiple (more than 2) marked locations. In the case of our problem, 
it is possible to reduce the multiple marked location case to the single location case at the cost 
of increasing the running time by a factor of logA'^ |AAn3j . Still, it would be interesting to be 
able to analyze the multiple item case directly. A recent paper by Szegedy jSzefl4j has shown how 
to analyze the multiple solution case for a different quantum walk algorithm, element distinctness 
|Ambn3j . It is open whether the methods from |Szen4j can be applied to search on grids. 
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